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Abstract 

The long-standing issue of the nature of the critical line of lattice QCD 
with the Wilson quark action at finite-temperatures, defined to be the line of 
vanishing pion screening mass, and its relation to the line of finite-temperature 
chiral transition is examined. Analytical and numerical evidence are presented 
that the critical line forms a cusp at a finite gauge coupling, and that the line 
of chiral transition runs past the tip of the cusp without touching the critical 
line. Implications on the continuum limit and the flavor dependence of chiral 
transition are discussed. 
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Elucidating the nature of chiral transition separating the high-temperature quark- 
gluon plasma phase from the low-temperature hadron phase has been one of the focal 
points of effort in recent numerical simulations of full lattice QCD including light 
dynamical quarks. Much of work in this direction has employed the Kogut-Susskind 
quark action since it retains a U(l) subgroup of chiral symmetry. On the other hand, 
studies with the Wilson quark action are less well developed in spite of the effort 
over the yearsjjl], |2|, |3|. In fact, two issues, which are fundamental for understanding 
the chiral transition with this action, have not been fully clarified to date. Both 
of the issues originate from the difficulty of identifying the chiral limit of massless 
quark due to the presence of explicit chiral symmetry breaking in the action, which 
is introduced to avoid the well-known species doubling problem. 

For the case of zero temperature pion mass has been conventionally employed 
to deal with this problem. Let (3 = Q/g 2 be the inverse gauge coupling and K 
the hopping parameter which controls quark mass. There is much evidence from 
analytical considerations and numerical simulations of hadron masses that the pion 
mass vanishes along a line K = K c (f3), the critical line, which runs from K c (/3 = 0) = 
1/4 in the strong-coupling limit to K c (/3 = oo) = 1/8 in the weak-coupling limit. At 
finite temperatures, one is then naturally led to define the critical line in terms of 
the pion screening mass extracted from the exponential decay of pion propagator for 
a large spatial separation. The first issue in understanding the chiral transition is 
how the critical line thus defined behaves as a function of (3 and the temporal lattice 
size N t , and how it is related to the critical line at zero temperature. 

The second issue concerns the relation of the critical line and the line of finite- 
temperature chiral transition K = K t (/3), the thermal line. Naively one would expect 
the thermal line to meet the critical line at some value (3 = (3 C . In this case, the 
region bounded by the thermal line and the critical line for j3 < j3 c represents the 
low-temperature phase and the other region the high-temperature phase. However, 
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initial simulations JT[ failed to find clear indications of such a behavior; the two lines 
appeared to run almost parallel toward the region of strong coupling, down to [3 = 
5.0 — 4.5 for the case of Nf = 2 flavors. These results raised the question if the 
two lines meet at all, which has led to subsequent studies[@, |J. In particular the 
QCDPAX Collaboration carried out simulations down to the strong coupling limit 
(3 = and concluded that the two lines meet at (3 ~ 3.9 — 4.00. However, their 
result has not satisfactorily answered the question in our view since the nature of the 
critical line at finite temperatures we have addressed above has not been clarified in 
their work. 

In this article we report results of our study on the two fundamental issues sum- 
marized above. Our analyses are based on the idea of spontaneously breakdown 
of parity and flavor as a characterization of the critical line, which has been put 
forward by one of usQ. The phase structure we found exhibits some unexpected 
features which we substantiate with analytical arguments and hybrid Monte Carlo 
simulations of Nf = 2 flavor full QCD with the Wilson quark action on an 8 3 x 4 
lattice. We shall also discuss implications of our results on the continuum limit of 
the Nf = 2 chiral transition and extensions for the case of Nf > 3. 

The starting point of our analysis is the Gross-Neveu model in two dimensions 
formulated with the Wilson quark action || . Except for confinement, this model 
has the features quite similar to those of QCD: asymptotic freedom, spontaneous 
breakdown of chiral symmetry and its restoration at a finite temperature. In the 
large N limit, the pion mass in this model is analytically calculable in the saddle- 
point approximation. In Fig. p] we plot the critical line corresponding to m n = on 
the (g, m) plane for the temporal lattice sizes N t = 2, 4, 8, 16 and oo where g and m 
are the bare coupling and quark mass, the latter related to the hopping parameter 
through K = I /(2m + 4). The result for N t = oo, which has been known for some 
time[|]], shows that the critical line forms three cusps which touch the weak-coupling 
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limit g = at m = 0, —2 and —4. The region outside of the critical line is a normal 
phase with m„ ^ 0, while the inside is the phase with spontaneous breakdown of 
parity characterized by < ^yy^ >^ 0. Conventionally one constructs the continuum 
limit of the model close to the point {g,m) = (0,0) along the branch of the critical 
line extending over oo > g > 0. The existence of additional branches converging 
toward (g,m) = (0, —2) and (0, —4) arises from a rearrangement of the doubler and 
massless fermion spectrum as m is varied. 

For understanding the phase structure at finite temperatures, the crucial feature 
revealed in Fig. [I] is that the three cusps of the critical line retract from the weak- 
coupling limit g = for finite temporal lattice sizes, forming a single continuous line 
which shifts toward strong coupling as N t becomes smaller. Thus, for a finite N t , 
the critical line is absent for sufficiently weak coupling. Also noteworthy is the fact 
that, even in the range of g where the critical line exists, the position of the critical 
line depends on Nt, albeit only slightly for large values of N t . 

The close resemblance of the Gross- Neveu model and QCD suggests that a similar 
behavior of critical line holds for QCD except that the number of cusps should 
increase to five because of the difference in the number of dimensions. For the zero 
temperature case, the available evidence supporting this possibility is as follows: 

(i) In the strong-coupling limit, analytical results in the large N limit for chiral 
observables such as < ip^Tzip > and associated susceptibilities show spontaneous 
breakdown of parity and flavor for \K\ > 1/4 @] . Results of a numerical simulation for 
the Nf = 2 full QCD system show good agreement with the analytic predictions 0. 

(ii) A quenched calculation of the number of conjugate gradient iterations needed for 
inverting the Wilson quark matrix over the range f3 = 5.0 — 6.0 revealed the presence 
of three peaks in the region K > for (3 > 5.60. (iii) A measurement of hadron 
masses at f3 = 6.0 in quenched QCD on a 10 3 x 20 lattice found clear evidence for 
three more values of K c beyond the conventional one at which pion mass vanishes 0. 



4 



The parity broken phase situated between two successive critical values is narrow, 
being of order 5K C « 0.05 — 0.1. 

For the finite temperature case, an apparent disappearance of the critical line 
toward weak coupling has been noticed in previous studies. In the Nf = 2 results of 
the MILC Collaboration H on an N t = 4 lattice in the range (3 = 4.9 — 5.3, the pion 
screening mass initially decreases toward larger values of K, but it increases beyond 
the thermal line K t ((3). A similar behavior is seen in the data of the QCDPAX 
Collaboration at /3 = 4.5 — 4.30. On the other hand, the conventional critical line 
has been shown to remain at (3 = and 3.50. Therefore, confirming a structure 
similar to that of Fig. |l]for QCD requires evidence that the conventional critical line 
turns back toward strong coupling at some value of f3. 

In order to examine this point, we have carried out hybrid Monte Carlo sim- 
ulations of Nf = 2 full QCD with the Wilson quark action on an 8 3 x 4 lattice. 
Runs were made in the range 3.0 < (3 < 5.3 and 0.15 < K < 0.31, which covers 
the region beyond the conventional critical line, typically in steps of 8(3 = 0.5 and 
5K = 0.01 — 0.02. The conjugate gradient method was employed for inverting the 
Wilson quark matrix with the stopping condition of ||£ — D*Dx\\ 2 /N < 10~ 5 with 
iV = ||x|| 2 or N = 12V with V the lattice volume. In order to maintain the ac- 
ceptance at the level of 70 — 80% or more, the hybrid Monte Carlo step size was 
decreased from 5t = 0.02 to 0.00125 as (3 is decreased or K is moved closer to the 
critical line. Choosing the interval of 0.5 or 0.25 to be the unit trajectory, at least 
50 — 100 trajectories were generated after thermalization with local observables mea- 
sured for each trajectory. At (3 — 4.0 and 3.5 hadron propagators were evaluated 
at every two trajectories by periodically doubling the lattice in spatial directions 
8 3 x 4 -> (8 x 2) x 8 2 x 4. 

In Fig. ^(a) we present our result for ir and p screening masses at (3 — 3.5 as a 
function of 1/K. Also plotted is the quark mass defined via the axial vector Ward 
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identity]?]]. To the right of the figure is the low-temperature phase bounded by the 
conventional critical line at K c w 0.2267 where the pion mass squared and the quark 
mass linearly vanish almost simultaneously. Clearly there exists another critical line 
at K c ~ 0.2454 below which the pion and quark masses behave in a reverse manner. 
Fig. [|(b) shows how the behavior changes at (3 — 4.0. We observe that the gap 
between the two critical values has either become extremely narrow or disappeared. 
These results lead us to believe that we have identified one of the cusps of the critical 
line, with the turning point located in the vicinity of /3 = 4.0. 

Let us now discuss the question whether the thermal line K t ((3) crosses the critical 
line K c {(3). In the light of our results above, it should be clear that this can not 
happen with our definition of K c (f3)\ the part of the critical line beyond the crossing 
point, assuming it exists, should be in the high temperature phase, where we expect 
the pion screening mass to be finite, while it vanishes along the entire critical line by 
definition. In other words the region close to the critical line has to be in the low- 
temperature phase even after it turns back toward strong coupling. This means that 
the thermal line should run past the turning point of the critical line and continue 
toward larger values of K. 

The validity of this consideration is confirmed through an examination of ther- 
modynamic observables. In Fig. |3|(a) we plot the real part of the Polyakov line and 
quark and gluon entropy densities in lattice units as a function of 1/K at (3 = 3.5. 
The two vertical lines show the position of the critical line estimated from the pion 
screening mass. For small values of 1 / K the observables take large values typical of 
the high temperature phase. However, toward the critical line at K c « 0.2454, they 
decrease, becoming roughly similar in magnitude to those on the other side of the 
conventional critical line at K c ~ 0.2267, which is in the low-temperature phase. 

Results for (3 = 4.0 are shown in Fig. 0(b). The vertical line marks the point 
where the two linear extrapolation of the pion screening mass squared in Fig. 0(b) 
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cross each other. The increase of the three quantities across the line shows that the 
thermal line comes close to the critical line at (3 = 4.0. 

In Fig. f| we summarize our findings for the phase structure of lattice QCD with 
two flavors of dynamical Wilson quarks on an N t = 4 lattice. Solid circles represent 
values of K c {(3) extracted from pion screening mass. We have tentatively assumed 
that the critical line turns back at /3 — 4.0. Open circles are estimates obtained by 
extrapolating the inverse of the number of conjugate gradient iterations in hybrid 
Monte Carlo runs to zero. The solid line, smoothly interpolating circles, represents 
one of the cusps of the finite-temperature critical line. (We expect the tip of the cusp 
to be actually rounded. See Fig. [l].) The thermal line is indicated by the dotted line, 
which is an interpolation of estimates from previous work (open squares) |], 0- @], 
continued toward larger values of K following our results. As we have argued the 
low-temperature phase goes around the tip of the cusp and extends into the region 
close to the upper part of the critical line. 

We emphasize that the zero-temperature critical line is absent from the point of 
view of finite-temperature partition function. For this reason we have not drawn 
the line in Fig. |j. We note, however, that available results indicate that this line 
practically forms a smooth continuation of the lower part of the finite-temperature 
critical line toward weak coupling. 

Let us comment here on the report of the QCDPAX Collaboration that the 
thermal line crosses the critical line at (3 ~ 3.9 — 4.0||. In our terminology their 
result is based on simulations carried out along the zero -temperature critical line 
from weak toward strong coupling. In view of the smooth connection of the zero- 
and finite-temperature critical lines noted above, the thermal line has to cross the 
zero-temperature critical line, which is the phenomenon suggested by the QCDPAX 
Collaboration. In this sense their result is consistent with our phase diagram. How- 
ever, we emphasize that this crossing point does not correspond to a singularity of 
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the partition function for a finite temporal lattice size. 

So far we have examined the case of N t = 4. For larger values of N t , we expect the 
thermal line and the cusp of the critical line to move toward larger f3 with the latter 
closely following the zero-temperature critical line. The distance between the thermal 
line and the tip of the cusp will diminish probably as 0(a) pa 0(1 /Nt). An important 
point to note is that the presence of explicit chiral symmetry breaking and simulation 
results including ours strongly indicate that the thermal line for the Nf = 2 case is a 
continuous crossover and not a true phase transition ||. Thus, a second-order chiral 
phase transition, as suggested by continuum sigma model analyses for this case|J], 
would emerge only in the continuum limit Nt — > oo. 

Another interesting question is how our results are modified for the case of Nf > 
3. Previous simulations support a first-order chiral phase transition || in agreement 
with sigma model analyses. This implies that the thermal line turns into a line of 
first-order phase transition near the tip of the cusp. The line will not cross the 
critical line for a finite temporal lattice size since our reasoning for the Nf = 2 case 
also applies for Nf > 3. 

Let us finally add that we expect the qualitative features of the phase structure we 
found to hold for a wider class of Wilson- type actions including the clover act ion [PHI. 
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Figure 1: Critical line for the lattice Gross-Neveu model on (g,m) plane. Temporal 
lattice size equals N t = 2, 4, 8, 16 and oo from inside to outside. 
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Figure 2: n and p screening masses and quark mass for Nf = 2 full QCD with 
Wilson quark action obtained on an 8 3 x 4 lattice periodically doubled in one of 
spatial directions. 
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Figure 3: Real part of Polyakov line(i?ef2), gluon (s g ) and quark (s q ) entropy density 
in lattice units for Nf = 2 full QCD with Wilson quark action obtained on an 8 3 x 4 
lattice. 
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Figure 4: Phase structure of Nf = 2 lattice QCD with Wilson quark action. Es- 
timates of the critical line for N t — 4 (solid line) and the thermal line for chiral 
transition (dotted line) are shown. For further details, see text. 
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